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Abstract 

In this work, we confront ten cosmological models with observational H(z) data. The possible in- 
teraction between dark energy and dust matter is allowed in some of these models. Also, we consider 
the possibility of (effective) equation-of-state parameter (EoS) crossing — 1. We find that the best 
models have an oscillating feature for both H(z) and EoS, with the EoS crossing —1 around redshift 
z ~ 1.5. 
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1 Introduction 

A lot of cosmological observations, such as SNe la [HE], WMAP [3], SDSS [2J, Chandra X-ray Observa- 
tory [5] etc., find that our universe is experiencing an accelerated expansion. These results also suggest 
that our universe is spatially flat, and consists of about 70% dark energy with negative pressure, 30% 
dust matter (cold dark matters plus baryons), and negligible radiation. Dark energy study has been one 
of the most active fields in modern cosmology [5]. 

In the past years, many cosmological models are constructed to interpret the present accelerated 
expansion. One of the important tasks is to confront them with observational data. The most frequent 
method to constrain the model parameters is fitting them to the luminosity distance 

which is an integral of Hubble parameter H = a/a, where a — (l+z) -1 is the scale factor (z is the redshift); 
a dot denotes the derivative with respect to cosmic time t. In this work, we use the observational H{z) data 
directly, rather than the luminosity distance di,{z). The observational H(z) data are based on differential 
ages of the oldest galaxies (7J. In [S], Jimenez et al. obtained an independent estimate for the Hubble 
constant by the method developed in [7J, and used it to constraint the equation-of-state parameter (EoS) 
of dark energy. The Hubble parameter depends on the differential age as a function of redshift z in the 
form 

H(z) = --±-%. (2) 
1 + z at 

Therefore, a determination of dz/dt directly measures H{z) [3]. By using the differential ages of passively 
evolving galaxies determined from the Gemini Deep Deep Survey (GDDS) [TD] and archival data [IT], 
Simon et al. determined H(z) in the range < z < 1.8 [5J. The observational H(z) data from [5] are given 
in Table Q] and shown in Figs. SHU 
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0.09 


0.17 


0.27 


0.40 


0.88 


1.30 


1.43 


1.53 


1.75 


H(z) (km s^Mpc- 1 ) 


69 


83 


70 


87 


117 


168 


177 


140 


202 


la uncertainty 


±12 


±8.3 


±14 


±17.4 


±23.4 


±13.4 


±14.2 


±14 


±40.4 



Table 1: The observational H(z) data [HE] (see [T5] also). 

These observational H(z) data have been used to constrain the dark energy potential and its redshift 
dependence by Simon et al. in [9] . Yi and Zhang used them to constrain the parameters of holographic dark 
energy model in |12j . Some relevant works also include [T2103]. Recently, in [T5], Samushia and Ratra have 
used these observational H(z) data to constrain the parameters of ACDM, XCDM and 0CDM models. 
In this work, we will revisit these observational H(z) data and compare them with some cosmological 
models. 

By looking carefully on the observational H{z) data given in Table [TJ and shown in Figs. [2H5J we 
notice that two data points near z ~ 1.5 and 0.3 are very special. They deviate from the main trend 
and dip sharply, especially the one near z ~ 1.5; the H(z) decreases and then increases around them. 
This hints that the effective EoS crossed —1 there. This possibility has not been discussed in previous 
works (e.g. [12], I15j). In the present work, we will seriously explore this possibility. 
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Figure 1: The quantity L(z) = H 2 (z)/H 2 — f2 m o(l + z) 3 versus redshift z, for the fiducial parameters 
H = 72 km s" 1 Mpc" 1 and fl m0 = 0.3. 



On the other hand, we also consider the possible interaction between dark energy and dust matter. 
This is inspired by the data point near z ~ 1.5, which dips so sharply and stays clearly outside of 
the best-fit of the ACDM, XCDM and 0CDM models studied in [TS]. In Fig. [Tj we show the quantity 
L(z) = H 2 (z)/Hq — + z ) 3 versus redshift z, which is associated with the fractional energy density 

of dark energy, for the fiducial parameters Ho — 72 km s _1 Mpc -1 and f2 TO o = 0.3, where the subscript 
"0" indicates the present value of the corresponding quantity. It is easy to see that the fractional energy 
density of dark energy of the point near z <~ 1.5 is negative (beyond ler significance). To avoid this, 
one can decrease the corresponding £l m o or make the matter decrease with the expansion of our universe 



slower than 



Inspired by this, it is natural to consider the possibility of exchanging energy between 



dark energy and dust matter through interaction. 

We assume that dark energy and dust matter exchange energy through interaction term C, namely 



pde + 3i? {pde + Pde) = ~C, 

p m + 3Hp m = C, 



which preserves the total energy conservation equation p to t + 3H (ptot + Ptot) 
extensively considered in the literature (see [rp ] fT7 1 [TH l fT ^ [2l ^ [2T | [2"2 l [2"3 l [2"i l l 

C (X Hp m , Hptot, Hpde, np m 4>, . . . 



(3) 
(4) 

0. The interaction forms 
1 [2S1 El] for instance) are 



In this work, we consider the simplest case for convenience, i.e. 

C = 3aHp m , 

where a is a dimensionless constant. Combining Eqs. (HI) and |(SJ), it is easy to get 

p m = p m0 a- 3 ^ = p m0 (l + z)^- a \ 



(5) 



(G) 



In the following sections, we will compare the observational H(z) data with some cosmological models 
with/without interaction between dark energy and dust matter. We consider a spatially flat FRW universe 
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throughout. We adopt the prior Hq = 72 km s _1 Mpc -1 , which is exactly the median value of the result 
from the Hubble Space Telescope (HST) key project [27] > an d is also well consistent with the one from 
WMAP 3-year result [3J. Since there are only 9 observational H(z) data points and their errors are fairly 
large, they cannot severely constrain model parameters alone. We perform a \ 2 analysis and compare the 
cosmological models to find out the one which catches the main features of the observational H (z) data. 
We determine the best-fit values for the model parameters by minimizing 

9 2 

X 2 (parameters) - ^ ^ Hmod ^ parameterS '' Z ^ ~ H °bs{zj)] ^ 
i=i 0-2 

where H m od is the predicted value for the Hubble parameter in the assumed model, H bs is the observed 
value, a is the corresponding lex uncertainty, and the summation is over the 9 observational H(z) data 
points at redshift Zi. 



2 ACDM and XCDM models without /with interaction 

In the spatially flat ACDM model, the Hubble parameter is given by 

H(z) = H v^Wl + zf + (1 - a m0 ) , (8) 

where V, m0 = K 2 p m o/(3-ffg) is the present fractional energy density of the dust matter, and k 2 = 8nG. 
By minimizing the corresponding x 2 , we find that the best-fit value for model parameter is fl m o = 0.302, 
while Xmin = 9-04 and Xmin/dof = 1-13, where dof is the degree of freedom. 

Then, we consider the interacting ACDM model (IntACDM). In this case, Eq. ([3]) becomes 

p A = -3aHp m . (9) 

By using Eq. ©, it is easy to find 

PA = Pmo a- 3(1 - a *> + const. (10) 

1 — a 

where const, is an integral constant. Inserting into the Friedmann equation H 2 = K 2 (p m + Pa)/3 and 
requiring H(z = 0) = one can determine the integral constant. Finally, in the IntACDM model, the 
Hubble parameter is given by 

H(Z) = H >h^(l + *)3(l-a) + ( 1 - p±) (11) 

y 1 - a \ 1 — a J 

By minimizing the corresponding x 2 , we find that the best-fit values for model parameters are £l m o = 0.386 
and a = 0.138, while X 2 mm = 8-89 and xLJdof = 1.27. 

We present the observational H(z) data with error bars, and the theoretical lines for ACDM (solid 
line) and IntACDM (dashed line) models with the corresponding best-fit parameters in the left panel of 
Fig.H 

In the spatially flat XCDM model, the Hubble parameter is 

H(z) = Eo/sUl + zf + (1 - n m0 )(l + z)3(i+«-x) , (12) 

where wx is the time- independent EoS of dark energy. By minimizing the corresponding % 2 , we find 
that the best-fit values for model parameters are fi m o — 0.284 and wx = —0.899, while Xmin — 9-^2 and 
xLn/ dof =1.29. 
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Figure 2: The observational H(z) data with error bars, and the theoretical lines for ACDM (left panel) 
and XCDM (right panel) models with the corresponding best-fit parameters for the cases without (blue 
solid line) and with (red dashed line) interaction, respectively. 



Next, we consider the interacting XCDM model (IntXCDM). In this case, Eq. ^ reads 

px + 377(1 + w x )px = -3aHp m . (13) 

Considering Eq. (J6j) , we obtain 

p x = const. a-^ 1+w ^ - a^- a \ (14) 

a + wx 

where const, is an integral constant. Again, inserting it into the Friedmann equation H 2 — K 2 (/0 m +Px)/3 
and requiring H(z = 0) = Ho, we can determine this integral constant. Finally, in the IntXCDM model, 
the Hubble parameter is written as 

H{ Z ) = H 0S [™X^°(1 + *)3(1-") + ( 1 - JH^hal) (1 + 3)3(1+**) . (15) 

w Y a + w x \ a + wxj 

By minimizing the corresponding \ 2 , we find that the best-fit values for model parameters are tt m o = 0.718, 
w x = -3.705 and a = 0.302, while xLin = 8 - 48 and X 2 mm l do f = lAl - 

We present the observational H(z) data with error bars, and the theoretical lines for XCDM (solid 
line) and IntXCDM (dashed line) models with the corresponding best-fit parameters in the right panel of 

Fig.m 

It can be seen clearly from Fig. [5] that none of the above four models may reproduce the sharp dip 
around z ~ 1.5; the data point near z ~ 1.5 deviates from model fitting by about 2a. 
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3 Vector- like dark energy 

As mentioned in the introduction, we are interested to seek a cosmological model whose effective EoS can 
cross —1. In this section, we consider the (interacting) vector-like dark energy model proposed in [28j (see 
also [TB]). In this model, the EoS of vector- like dark energy and then the effective EoS can cross —1 in 
principle. From [161 [28] , the energy density and pressure of vector-like dark energy are given by 



3 

PA = ~ 2 

PA = 



(a + HA\ + W(A 2 ), (16) 
\{a + HA) 2 -W(A 2 )+2^A 2 , (17) 



where A 2 (t) is the time-dependent length of the so-called "cosmic triad" of three mutually orthogonal 
vector fields. In this work, we would like to consider the case with exponential potential, namely 

V (A 2 ) = V A exp (-Xk 2 A 2 ) , (18) 

where Va and A are constants. By using Eq. (0) and new quantities B = kA and Va = k 2 Va, we can 
recast Eq. §S§ as 

(b + HB^ B + 'iHB + [h + 2H 2 ^j B -2\V A Be- XB2 = -3aHHfa m0 (l + z) 3(1 " q) . (19) 
To find the solution for B(t), we can use the initial conditions 



/2(A-1) 

Bini — ±Y -r , B ini — 0, (20) 

which come from |16j . under the conditions A > 1 and a < 1. It is convenient to change the time t to 
redshift z, and 

f=-(l + z)Hf, (21) 
/ = (1 + z)H 2 f - (1 + z)Hf + (1 + z) 2 H 2 f" , (22) 

where /' = df/dz for any function /; we have used Eq. (J2J) . From the Raychaudhuri equation H = 
-k 2 (pa + p m +pa)/2, we get the H in Eq. ([22| as 

H=- [-(1 + z)B' + B} 2 H 2 + XV A B 2 e- XB2 - ^H 2 fl m0 (l + z) 3 ^. (23) 
From the Friedmann equation, the Hubble parameter reads 



H = H Jn m0 (i + z ) 3 ( 1 -°>) + ,, " /M 



Inserting Eq. (116p into it, we find 



3H 2 



H 2 = H 2 n m0 (l + z) 3{1 - a) +VAe XB 
\-\[-{l + z)B' + B] 2 

From Eq. ([M]) and the requirement of H(z = 0) = Hq, only three of a, A, TO o and Va are independent 
of each other. We choose our free model parameters as a, A and f2 TO o. 

Now, we can numerically solve Eqs. (JTSJ) and with the help of Eqs. ([2T 1) -([M] 1 to get the B(z). 
Then, we obtain the Hubble parameter H(z) from Eq. (f2"ll) . So, the corresponding \ 2 is in hand. 



7 




0.25 0.5 0.75 1 1.25 1.5 1.75 0.25 0.5 0.75 1 1.25 1.5 1.75 

z z 

Figure 3: The left panel is the observational H (z) data with error bars, and the theoretical lines for VecDE 
(blue solid line) and Int VecDE (red dashed line) models with the corresponding best-fit parameters. The 
right panel shows the effective EoS (thick lines) and the EoS of vector-like dark energy (thin lines), for 
VecDE (blue solid line) and Int VecDE (red dashed line) models. 



For the vectordike dark energy model without interaction (VecDE), i.e. a = exactly, we find the 
minimal \ 2 as xLn = 9-05 (Xmin/dof = 1-29), for the best-fit parameters A = 1.001 and f2 m o = 0.299 
(the corresponding Va = 0.70 Hq). 

For the interacting vector-like dark energy model (Int VecDE), i.e. leaving a as a free parameter, we 
find the minimal % 2 as Xmin = 9-04 (Xmin/dof = 1-51), for the best-fit parameters a — —0.002, A = 1.005 
and f2 m o = 0.30 (the corresponding Va = 0.703 Hq). 

We present the observational H(z) data with error bars, and the theoretical lines for VecDE and 
Int VecDE models with the corresponding best-fit parameters in the left panel of Fig. [3] The right panel 
of Fig. [3] shows their corresponding effective EoS and the EoS of vector-like dark energy. 

We see that the H(z) predicted by VecDE and IntVecDE models, and their corresponding effective 
EoS, cannot be clearly distinguished. Also, it can be seen clearly from Fig. [3] that none of these two 
models may reproduce the sharp dip around z ~ 1.5; the data point near z ~ 1.5 deviates from model 
fitting by about 2a. Although the EoS of vector-like dark energy for VecDE and IntVecDE models are 
different, they have not crossed —1. We will briefly discuss this point in section \5\ 

4 The models with oscillating H(z) ansatz 

Obviously, all six models studied above fail to catch the features of the sharp dip around z ~ 1.5 and EoS 
crossing —1 mentioned in the introduction. In this section, we consider some parameterized models. 
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The first parameterized model (OA1) is the best model studied in [SHj which fits the SNe la data very 
well. Its Hubble parameter is given by 



H(z) = H y/Sl m0 (l + z) 3 + ax cos(a 2 z 2 + a 3 ) + (1 - ai cosa 3 - f2 m0 ) , (25) 

where a\, a 2 and a-$ are constants. By minimizing the corresponding \ 2 for the H{z) data, we find that 
the best-fit values for model parameters are fl m o = 0.241, oi = 1.316, a 2 = 2.717 and a$ = —3.933, while 
xLn = 4-27 and xLJdof = 0.85. 

Adding interaction into it, we obtain the IntOAl model with Hubble parameter 



H{z) = H ^Q m0 (l + z) 3 ( 1 -«) + a x cos(a 2 2 2 + a 3 ) + (1 - ai cosa 3 - Q m0 ) , (26) 

We find that the best-fit values for model parameters are £l m o — 1-002, a — 0.407, a\ = —1.419, a 2 = 2.995 
and 03 = —1.366, while Xmin ~ 3-54: and Xmin/dof = 0.89. 

We present the observational H(z) data with error bars, and the theoretical lines for OA1 and IntOAl 
models with the corresponding best-fit parameters in the left panel of Fig. [U The right panel of Fig. 0] 
shows the effective EoS, 

WaffS B«=-l+hl + z)% 1 (27) 
Plot o ti 

for OAl and IntOAl models, respectively. It is obvious that both effective EoS of OA1 and IntOAl models 
crossed —1. The H(z) lines predicted by these two models dip around z ~ 1.5. So, it is not surprising 
that they fit the observational H(z) data much better than the six models studied above. 




Figure 4: The left panel is the observational H{z) data with error bars, and the theoretical lines for OAl 
(blue solid line) and IntOAl (red dashed line) models with the corresponding best-fit parameters. The 
right panel shows their corresponding effective EoS. 



9 



Next, we consider a variant (0A2) of the best model studied in [30] which also fits the SNe la data 
very well. Its Hubble parameter reads 

H(z) = H yJn mQ (l + zf + ai(l + zf [cos(a 2 2 2 + a 3 ) 1 cosa 3 ] + (1 - Q m0 ) , (28) 

where a%, a 2 and a 3 are constants. By minimizing the corresponding x 2 , we find that the best-fit values 
for model parameters are f} m o = 0.287, a\ — 0.132, a 2 = 2.971 and a 3 = —4.481, while Xmin ~ 2-81 and 
XmiJ do f = 0.56. 

We get the IntOA2 model by adding interaction into OA2. The corresponding Hubble parameter is 
given by 

H(z) = H a y/ fi m0 (l + z)3(!-«) + + z)3(i-°0 [cos(a 2 z 2 + a 3 ) - cosa 3 ] + (1 - fi m0 ) • (29) 

We find that the best-fit values for model parameters are tt m o = 0.252, a = —0.048, a\ = —0.117, 
a 2 = 2.975 and a 3 = -1.351, while xLin = 2 - 80 and Xmin/ do f = °- 70 - 

We present the observational H(z) data with error bars, and the theoretical lines for OAl and IntOAl 
models with the corresponding best-fit parameters in the left panel of Fig. [5] The right panel of Fig. [S] 
shows their corresponding effective EoS. Again, we see the OA2 and IntOA2 models catch the main 
features of the observational H(z) data; they reproduce the sharp dip around z ~ 1.5 and their effective 
EoS crossed —1 also. Since the value of a is too small, there is no significant difference between OA2 and 
IntOA2 models, unlike the case of OAl and IntOAl models. 




0.25 0.5 0.75 1 1.25 1.5 1.75 0.25 0.5 0.75 1 1.25 1.5 1.75 

z z 



Figure 5: The left panel is the observational H(z) data with error bars, and the theoretical lines for OA2 
(blue solid line) and IntOA2 (red dashed line) models with the corresponding best-fit parameters. The 
right panel shows their corresponding effective EoS. 
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5 Concluding remarks 

In Table [21 we summarize all ten models considered in this work. It is obvious that they divide into two 
groups. The first six models fail to catch the main features of the observational H(z) data, and hence 
have large Xmin anc ^ Xmin/d°f- The last four models with oscillating H(z) ansatz, on the contrary, are 
well compatible with the observational H(z) data. They can reproduce the sharp dip around z ~ 1.5, 
and their effective EoS crossed —1 also. Therefore, it is not surprising that their corresponding Xmin an d 
Xmin/ do f are significantly smaller than the ones of the first six models. 



Model 


Xmin 


X mm / d °f 


^ (x ^* Xmin) 


Ranked by 


ACDM 


9.04 


1.13 


0.34 


5 


IntACDM 


8.89 


1.27 


0.26 


6 


XCDM 


9.02 


1.29 


0.25 


7 


IntXCDM 


8.48 


1.41 


0.21 


9 


VecDE 


9.05 


1.29 


0.25 


8 


IntVecDE 


9.04 


1.51 


0.17 


10 


OA1 


4.27 


0.85 


0.51 


3 


IntOAl 


3.54 


0.89 


0.47 


4 


OA2 


2.81 


0.56 


0.73 


1 


IntOA2 


2.80 


0.70 


0.59 


2 



Table 2: Summarizing all ten models considered in this work. 

Some remarks are in order. The first one is on the failure of the vector-like dark energy model to 
reproduce the main features of the observational H{z) data, since the EoS of vector-like dark energy 
and then the effective EoS can cross —1 in principle. According to T6J, for the case of interaction term 
C = 3aHp m , the vector-like dark energy behaves like a pure cosmological constant in the late time. So, it 
is not surprising that the vector-like dark energy model considered in this work is very close to the ACDM 
model. However, when the interaction term C takes other forms, the late time behavior of vector-like dark 
energy can be considerablely different, as shown explicitly in |16j . Therefore, it is interesting to compare 
the observational H(z) data with the interacting vector-like dark energy model with other interaction 
forms and different potentials. We leave this to future works. 

Secondly, we note that the best-fit value fl m o = 1-002 of IntOAl model is inconsistent with the results 
from clusters of galaxies [3S] and 3-year WMAP [3] etc. In this sense, the IntOAl model may be ruled 
out, although it fits the observational H(z) data fairly well. 

The third remark is concerned with EoS crossing — 1. There are many pieces of other observational 
evidence for this in the literature [5U1 13H 1321 133] ■ Also, a lot of theoretical models whose EoS can cross 
— 1 have been built (see for examples [TIB [TTJ [2U [321 [3H HO] and references therein). In this work, we 
present independent evidence for EoS crossed —1, from the observational H{z) data. 

Fourthly, we see that the last four models fit the data much better than the first six models. We note 
that the last four models have an oscillating feature, for both predicted H(z) and EoS, as shown in Figs. 3] 
and [5) In fact, this is also the main ideas of [33 [3S1 133 IS] • However, the last four models studied here 
and the ones considered in [33 [3H [33 HJJ are all parameterized models. We consider that it is important 
to build some physically motivated models in which the oscillating EoS arises naturally. 
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Fifthly, we admit that although the first six models have considerably larger Xmin/dof, they are not 
unacceptable. For instance, the ACDM model has Xmin = 9-04 for eight degrees of freedom with 34% 
probability, which is not unacceptably low. Also, before the new and improved H(z) data are available, it 
is too early to talk about strong observational evidence for non-monotonic behavior of H(z), since one can 
wonder that the dip around z ~ 1.5 could be due to unknown measurement errors and so on. Therefore, 
to firmly rule out non-oscillating models (e.g. the first six models considered here), the more and better 
observational H(z) data are required. 

Finally, the observational H{z) data provide an independent approach to constrain the cosmological 
models. However, by now, the observational H(z) data only have 9 data points and their error bars are 
fairly large. Hence, they cannot severely constrain the cosmological models so far. A good news from [IS] 
is that a large amount of H(z) data is expected to become available in the next few years. These include 
data from the AGN and Galaxy Survey (AGES) and the Atacama Cosmology Telescope (ACT), and by 
2009 an order of magnitude increase in H(z) data is anticipated. 
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